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This  Part  1  of  a  projected  2-part  study  initiates  an  indepth  network  study 
of  the  inverse  problem  for  1-dimensional  passive,  lossless  smooth  disper¬ 
sive  nonuniform  structures.  The  diagnostic  procedure  is  via  time-domain 
pulse-reflectometry  and  rests  on  a  new  continucjus-discrete  inversion 
formula.  .  — 1 — 

It  is  expected  that  this  formula,  when  used  together  with  the  weak  conver¬ 
gence  principle  described  in  Section  V.  will  lead  to  a  well-posed  (over) 
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I.  INTRODUCTION 


It  has  only  recently  been  realized  in  engineering  circles,  largely  due 
to  the  outstanding  pioneering  work  of  Gopinath  and  Sondhi  [5,  6],  that  the 
solution  of  the  inverse  problem  for  the  classical  lossless  nonuniform  line 
possesses  a  simple  and  natural  time-domain  setting.  This  simplicity  is 
directly  attributable  to  the  fact  that  the  expression  of  causality  in  time  is 
incomparably  less  complicated  than  its  expression  in  terms  of  frequency 
behavior. 

Nevertheless,  the  problem  of  diagnostics  for  dispersive  passive, 
lossless,  smooth  1 -dimensional  tapers  remains  open  despite  the  fact  that 
its  discrete  counterpart  has  already  been  solved  [2,3]!  In  this  part  I  of  a 
projected  2 -part  study  we  demonstrate  that  the  results  obtained  in  Ref,  3 
yield  the  deep  continuous -discrete  inversion  formula  (47)  which  provides 
a  solid  network  foundation  for  our  entire  approach.  La  fact,  when  coupled 
with  the  weak  convergence  principle  described  in  Section  V,  it  immediately 
leads  to  an  effective  integral  equation  for  the  logarithmic  taper  P(x)  of  the 
telegrapher's  line  (theorem  2)  in  terms  of  the  input  impulse  response. 

We  should  like  to  stress  that  our  conceptualization  of  a  smooth  taper 
is  based  squarely  on  physical  considerations  and  the  role  of  dispersion  is 
given  a  quantitative  characterization. 

It  is  also  expected  that  the  final  solution  will  involve  a  pair  of  coupled 
Fredholm  integral  equations  of  the  second  kind.  Thus,  the  inverse  problem, 
in  its  full  generality,  appears  to  possess  the  essential  property  of  well- 
posedness.  Moreover-,  the  corresponding  numerical  algorithms  are  of  the 
Levinson  type  and  the  theory  of  orthogonal  polynomials  on  the  unit  circle 
can  be  exploited  to  advantage. 

As  regards  notation,  most  of  it  is  self-explanatory.  However,  we 

point  out  that  for  a  matrix  A,  A,  A7,  A*(sA7)  and  det  A  denote  its  complex 

conjugate,  transpose,  adjoint  and  determinant,  respectively.  As  usual, 

column-vectors  are  written  as  a,  b  or  as  x  =  (x, ,  x_, . . .  ,x  )7  if  it  is 

desirable  to  exhibit  its  components,  and  0,1,0  stand  for  the 

-n  n  m  9  n 

n-dimensional  zero  vector,  the  nxn  identity  matrix  and  the  mxn  zero 
matrix,  respectively.  We  also  write  A>0fi  to  indicate  that  the  hermitean 
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matrix  A=  A  is  nonnegative-definite.  Lastly,  if  A(p)  is  any  real  meromorphic 
matrix  function  of  p,  it  is  convenient  for  typographical  reasons  to  let 

\(p)  s  A/<-P>  • 

Clearly,  if  p  =  jU),  <s  real,  A^juu)  =  A*(ju)). 

Since  this  is  a  research  report  our  demands  on  the  reader's  knowledge 
vary  widely.  Nevertheless,  a  solid  exposure  to  the  positive-real  coocept 
and  its  applications  to  insertion-loss  design  [l  ]  should  suffice  for  the  most 
part. 


II.  THE  MODEL 


From  an  engineering  point  of  view  it  appears  reasonable  to  con¬ 
ceptualize  a  continuous  1 -dimensional  electrical  taper  (Fig.  la)  as  the 
"physical"  limit  of  a  discrete  cascade  of  generic  type  shown  in  Fig.  lb. 

If  the  taper  is  also  internally  passive  and  lossless  and  operates  in  a  single¬ 
mode  environment,  it  is  then  natural  to  impose  the  following  constraints. 

C^)  Every  one  of  the  n  LLN's  is  a  lumped  passive,  lossless  2-port. 

C^)  The  n+1  lines  are  all  ideal  TEM  and  possess  the  same  1-way 
delay  T>0  but  arbitrary  positive  characteristic  impedances  Rq,R^,  . . .  ,  R^. 
In  our  opinion,  the  presence  of  the  initial  line  (of  characteristic  impedance 
RQ)  is  indispensable  and  serves  to  provide  the  necessary  single-mode 
transition  between  the  actual  input  and  the  actual  start  of  the  taper.  (In  any 
case,  only  experiment  can  decide  for  or  against  this  assumption.) 

C^)  The  total  1-way  delay 

T  =  (n+l)T  (1) 


remains  constant  as  n  •»«. 


Let 


zm(p)  denote  the  impedance  seen  looking  into  the  input  side  of 


LLN  m  with  its  output  closed  on  R^  (Fig.  2a).  By  definition, 


i  \  _  Zm^  "  Rm-1 

sm(p)  ’  z(p)  +  Rm  ,  * 

m  m- 1 


m=l-»  n 


(2) 


is  the  corresponding  normalized  "junction"  reflection  coefficient.  Given 
that  the  load  termination  is  passive  but  unknown,  the  purpose  of  internal 
diagnostics  for  the  discrete  cascade  is  to  determine  as  much  as  possible 
about  these  coefficients  by  means  of  measurements  carried  out  exclusively 
from  the  input  side  of  the  structure. 

Now  if  the  notion  of  a  physical  limit  is  meaningful,  LLN  m  and  the 

line  of  characteristic  impedance  Rm  contiguous  on  its  right  must  constitute 

an  accurate  representation  of  the  properties  of  a  section  of  continuous  taper 

2 

situated  between  x=mr  and  x=  (m+l)T,  at  least  in  the  limit  of  small  t.  We 


*p=a+j(t)is  the  complex-frequency  variable. 

2 

The  analogy  with  the  discrete  model  requires  that  x  should  equal  the 

time  in  secs,  for  an  electromagnetic  pulse  injected  at  x=0  to  reach  the 
cross-section  identified  by  the  abscissa  x. 


3 


are  therefore  led  to  introduce  the  following  interpretive  assumptions. 

Aj)  There  exists  a  positive  and  continuously  differentiable  function 
R(x)  which  plays  the  role  of  a  “local"  characteristic  impedance  for  the 
continuous  taper. 

A^)  Let  z(x,  p;T)  denote  the  impedance  seen  looking  into  a  section 

of  continuous  taper  lying  between  x  and  x-fT  which  at  x+T  is  closed  on  the 

resistance  R(x+t)  (Fig.  2b).  Then,  there  exists  a  function  r|(x, p)  which 

is  continuously  differentiable  in  x  for  every  fixed  p  in  Rep>  0,  and  real 

3 

and  rational  in  p  for  every  fixed  x  such  that 

=  l+2TTi(x,p)  +  o(x,p;T)  ,  (3) 

0<x<  T  . 

A^)  For  mT-*x>0,  x  fixed,  and  for  p  fixed  in  Rep>0, 


z(x,  p;t)  -R(x)  _ 
z(x,  p;T)  +  R(x)  " 


sm(P>  +  o(x,p;T)  . 


(4) 


As  an  immediate  consequence  of  the  above  three  assumptions  it  is 
easily  shown  that  as  mT-»x>0,  x  fixed, 

s  (p) 

limit  —  =  P(x)  +  T|(x,  p)  (5) 

t|o  t 

where 


P(x) 


(6) 


is  the  so-called  logarithmic  taper.  The  two  functions  P(x)  and  T|(x,  p)  de¬ 
limit  the  generic  taper  profile  while  formula  (5)  supplies  the  means  for  their 
calculation  through  the  discrete  approximation. 

Our  first  objective  is  to  derive  an  expression  for  sm(P)  *n  terms  of 
input  data.  Since  interaction  with  the  load  must  be  avoided  it  is  natural 
to  seek  to  develop  a  procedure  based  on  time-domain  pulse -reflectometry. 


limit  P^xtf  =  0  ,  uniformly  in  x. 

rjo 


-4- 


HI.  DIAGNOSTICS  FOR  THE  DISCRETE  CASCADE 

Let  the  discrete  cascade  be  excited  through  the  input  current  driver 
i(t) :  6{t)  and  let  h^(t)  denote  the  corresponding  input  voltage  impulse 
response  v(t)  under  initially  quiescent  conditions. 

Evidently,  because  the  round-trip  travel  time  for  each  line  is  2T  secs., 

h^(t)  is  independent  of  the  load  for  all  t  in  the  interval  0<  t<  2(n+l)T.  In 

/n\  “ 

addition,  in  this  same  time  interval  h'  '(t)  admits  a  decomposition  into  a 
sum  of  "reflected"  waves  of  the  form, 


.  (n).  . 

=  6(t)  +  Yj  2hm(t-2mT) 
°  m=  1 

whose  properties  are  summarized  in  theorem  1. 
Theorem  1  (Appendix).  Let 

00 

„,..}<»)  =  J 

m  nn  *  _  m 


h  (t)«  c  (juj)  =  f  e”^U)t  h  (t)  dt  , 
m  '  mJ  '  J  _  m1 

-00 


(7) 


(8) 


and  let 


Am(p)  =  det 


m  =  l-*n 


1 

C1*<P> 

•  •  •  ctn*(p) 

c^p) 

1 

cl*(p) 

*  *  *  cm-l*^ 

c2(p) 

c^P) 

1 

•  *  *  Cm-2*^ 

• 

• 

• 

•  •  «  • 

• 

c  (p) 
m'r' 

• 

‘m-lW 

• 

•  •  *  • 

.  .  .  1 

m  =  0  -»n 


.  (9) 


Then, 


1)  every  Laplace  transform  cm(p)  is  a  real  rational  function  of  p 
analytic  in  Rep>0,  p=oo  included.  (Thus,  hm(t)  vanishes  for  t<  0,  m=  l-*n.) 

2)  With  the  possible  exception  of  a  finite  number  of  id's  at  most, 

A  (ju>)  >  0  ,  m  =  0  -»n  ,  (10) 

m 
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uu  real. 


3)  Let 


S(m)(p)  = 


1  9r»(P>  8\T^> 

■iT’w 


denote  the  scattering  matrix  o£  LLN  m  normalized  to  R  .  on  the  left  and 

/m\  m-x 

Rm  on  the  right,  let  dm(p)  =  det  Sv  '(p)  and  let  Am/  (p)  denote  the  minor  of 
Am(p)  that  results  by  suppressing  its  first  row  and  last  column.  Then,  ^ 

S1™)(P)  sm(P> 

—  =  -wr  *  (12) 

S?l*tP)  S7?*(P) 


3  (P)  = 


Am'<P>  T1  ,-l 


Am-1^  r=0 


.  n  d/(P)  , 

_ A  * 


Am(P>am-Z<P> 


m  =  l-»n;  dQ(p)  =A_L(p)  a  1  .  (15) 

To  help  the  reader  develop  an  appreciation  for  the  key  diagnostic 
formula  (13),  we  shall  use  it  to  solve  a  problem  of  considerable  interest 
in  its  own  right. 

Problem  1.  Determine  the  dielectric  constants  of  a  cascade  of  n+1 
nonmagnetic  ideal  TEM  lossless  layers  of  the  same  1-way  delay  T  by  means 
of  input-side  time-domain  reflectometry. 

Solution.  Clearly,  layer  m  acts  as  a  TEM  line  of  characteristic 
impedance 

R  =  ohms  ,  (16) 

m  sr 

m 

where  em  >  0  is  its  relative  dielectric  permittivity,  m  =  0-»n.  Further,  all 
LLN'  in  Fig.  '  collapse  into  a  pair  of  perfectly  conducting  wires  so  that 


-<f(p>  and  sj“)(p)  cannot  both  be  identically  zero  since  S^mHp)  is 
regular  pataunitary. 


S<m>(p>  - 


.  R  -R  ,  2  vmt  71 

1  m  m-1  tn  m- 1 

+  ftm-l  2vft  ft  ;  R  .  -R 

L  m  m-1  m-1  m 


and  dm(p)  =  -1,  m  =  l-»n.  Consequently,  the  c^'s  are  real  constants  and 


,  (n),  . 

ii^-^  =  6(t)+  V  2cmft(t-2m')  . 


Of  course,  R^  is  obtained  as  the  strength  of  the  initial  voltage  impulse  which 
appears  at  the  input  interminals. 

Hence,  in  this  all-transmission  line  case,  z  (p)  =  R  , 

’  *  m  r  m 

R  -R 
m  m-1 

sm  =  R— TE - -  <‘9> 

m  m-1 


R  =  R  •  -  ,  m  =  1  -» n  .  (20) 

m  m-1  1  -  s  ’  •  \  i 

m 

/  r»\ 

According  to  (18),  h'  (t)  is  the  sum  of  an  initial  impulse  of  strength 
Rq  and  n  successive  reflected  impulses  of  strengths  2RQcm  ^  separated  by 
2t  secs.  Thus,  all  cm's  can  be  measured  and  (13)  yields  the  explicit  answer, 

l  ^  ' 

sm  =  (’1}  •  £2L_  ’  m=  1_*n  ’  (21) 

m-1 

where 


C1  1 

A  =  det  c_  c. 

m  4  1 


1_C  m  cm-l  Cm-2  *  *  J 

and  is  the  minor  formed  by  striking  out  its  first  row  and  last  column. 
Formulas  (20)  and  (21)  constitute  the  desired  solution. 


For  t>  0  the  input  is  open-circuited  and  this  accounts  for  the  doubling, 


Tmzzx—Tza 


Before  proceeding  to  the  continuous  taper  it  is  necessary  to  clarify 
two  conceptual  points  of  great  importance. 

Point  1.  First,  we  must  show  that  the  coefficients  c__  (p),  m=  l-*n,  are 

"  i  i  1*1 

actually  determined  uniquely  by  the  segment  of  h  n  (t)  over  0<t<2(n+l)T.  To 
this  end,  observe  from  (7)  that 


h^(t) 

2h.  (t-2  T )  ,  2t  <t<  4t 
R0  1  — 


(23) 


and  hj(t)  is  therefore  known  for  0<  t<  2t.  Since  h^(t)  is  expressible  as  the 
sum  of  an  impulse  at  t  =  0  and  a  finite  number  of  polynomials  in  t  multiplied 
by  exponentials,  its  values  over  2T<t<  4t  fix  it  completely  for  all  t>0  by 
analytic  continuation. 


Similarly,  for  4t  <  t  <  6t  , 


2h2(t-4T)  =  iK^a.-2h1(t-2T) 


(24) 


which  means  that  h^(t)  is  also  known  for  0<  t<  2t  .  Hence,  by  the  above 
reasoning,  it  too  is  defined  uniquely  for  all  t>0,  etcetera. 

Point  2.  Second,  although  the  coefficients  c__(p),  m  =  l-»n,  are  unique 
constructs  from  the  given  data  h'  *(t),  0<t<  2(n+l)T,  the  application  of 
formula  (13)  also  requires  knowledge  of  the  n-1  functions  d^(p),  d^(p),  . . . , d^_^(p). 
In  a  great  number  of  interesting  situations  it  is  possible  to  ascertain  the  d's 
in  terms  of  Rq  and  the  c's  with  the  aid  of  the  following  simple  rule. 

Rule  (Appendix).  Suppose  that  each  LLN  is  all-pass  free  on  its  output 

“  £  1 

side  and  D.C.  transparent.  Let  the  strict  Hurwitz  polynomial  a^(p)  be  con¬ 
structed  with  all  the  zeros  of  ^(p)/A^(p)  in  Rep<0,  multiplicities  in¬ 
cluded.  Then, 

k  -i  k  Vp) 

“0d>-(-‘)  -w  ’  1251 


k  =  0-+  n- 1 ;  an(p)  s  1 


For  example,  if  k=l,  the  polynomial  a^(p)  is  constructed  with  the  closed 
left-half-plane  zeros  of 

LLN  m  is  D.C,  transparent  if  S^m^(0)  is  given  by  (17). 


-8- 


(27) 


-1  ftl(p) 

di  {P)  =  "  o^TpT  *  etc* 

Fortunately,  physical  tapers  are  usually  constructed  with  material 
that  cannot  produce  zeros  of  transmission  in  Rep>  0  and  their  all-pass  free 
character  is  intrinsic. 


IV.  DIAGNOSTICS  FOR  THE  CONTINUOUS  TAPER 


:(“)/ 


We  lay  down  the  premise  that  the  normalized  scattering  matrix  S'  '(p)  of  a 
section  of  continuous  taper  lying  between  x  =  mT  and  x=(m+l)r  approximates 
the  identity  transformation  as  T  -*  0.  That  is,  for  every  fixed  p  in  Rep>  0, 


’0 

r 

1 

0 

limit  S(m)(p)  = 
T^O 


Hence,  to  first  order  in  T  and  for  Rep>  0, 

A^m)(p)  =  l2-TZ(x,p) 

where 


‘  0  1 
1  0 

■ 


and 


Z(x, 


’  zn(x.p) 

Zi2(x,p) 

z21(*.p) 

*22(x.p) 

is  a  real  rational  matrix  in  p  for  every  fixed  x  in  0<  x<  T. 

Since  the  taper  is  passive, 

l2-S(m)*(p)S(m)(p)  =  T  (Z(x,  p)  +  Z*(x,  p))  >  02 

for  all  p  in  Re  p>  0.  It  follows  that 

Z(x,  p)  +  Z^x,  p)  >  02  , 

0  <  x  <  T  ;  Re  p>  0  , 

Again,  because  the  taper  is  also  lossless, 

l2-s£m)(p)S(m)(p)  =  T(Z(x,p)  +  Z#(x,p))  «  02 

and  Z(x,  p)  is  forced  to  be  skew- parahermitean  in  p: 


(28) 


(29) 


(30) 


(31) 


(32) 


(33) 


(34) 


The  2x2  matrix  Z(x,  p)  supplies  the  taper's  per-unit  length  description. 
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Z(x,p)  +  Z„,(x,p)  s  02  ,  (35) 

0<  x<  T;  all  p 

In  short,  for  the  class  of  passive  lossless  continuous  tapers  under  considera¬ 
tion,  the  per-unit  2x2  descriptor  Z(x,  p)  is  a  real  rational  Foster  matrix  in 
p  for  every  fixed  x.  In  particular,  zn(x,p)  and  z22(x,p)  are  reactances  and 

z21(x.  p)  =  ‘Z12*(*.P)  •  (36) 

From  Eqs.  (29)  and  (31), 


s(m)(P) 


*-tz21(x,  p)  1-tz22(x,; 

l-Tzll(x.P)  -Tzu(x,] 


Hence,  a  comparison  with  (11)  yields  the  first-order  identifications, 

sm(p)  =  “Tz21(x'p)  * 

S^^p)  =  1  -  ^z22(x,p)  =  exp[-Tz22(x,p)]  , 

4T)<p> =  1  -tzu(x»p)  =  exp|^- Tz j  j (x,  p) J  , 

S(2™l(p)=  -TZn(x,p)=  -  Sm4,(p) 


d  (p)  = 


=  -exp[-2T§(x,  p)] 


where 


?(x,p)  = 


zn(x,p)  +  z22(x,p) 


Clearly,  §(x,  p)  is  a  real  rational  reactance  function  of  p  for  every 
fixed  x.  Moreover,  as  is  readily  appreciated  from  (42), 


If  the  taper  is  reciprocal,  s'™ '(p)  = 
and  ?(x,p)  =  zn(x,p). 


s2™;(p)  so  that  zn(x,p)  =  z22(x,  p) 


\ 


(44) 


Tg(x,OJ) 


dp 


p=j(l) 


represents  the  average  to-and-fro  group  delay  per  second  of  taper  at  location 
x  due  to  dispersion. 

We  can  now  recognize  with  the  aid  of  (42)  that  for  tJo  and  x  =  mT  fixed, 


m-1 


m-1 


n  d"1(p)  =  (-l)R1"1  .  exp^T  J  ?(rT,p)]+  O(T) 
r=0  r=0 


(45) 


=  (-l)1”'1  •  exp^2  J  ?(y,  p)dy]  +  O(t) 


(46) 


Consequently,  since  s  (p)/r  =  -z7.(x,p),  Eqs,  (15),  (13)  and  (46)  yield  the 

m  m  Jt 

fundamental  c ontinuous -di  s  c  r ete  inversion  formula. 


P(x)  +  T)(x,p) 


'21 


(x,p)  =  exp[2  j  ?(y,p)dy] 


limit  (-l)m~* 

t|0 


A  /(p) 
m 


‘  TAm-l(P'  ' 


x=mT  fixed;  0<x<T  .  (47) 

In  general,  z^fx,  p)  determines  only  the  sum  P(x)+T|(x,  p)  and  not  P(x) 
and  ^(x,  p)  individually.  However,  if  the  taper  is  D.  C.  transparent  so  that 
it  reduces  to  a  pair  of  perfectly  conducting  wires  at  uu=0,  then  T|(x,  0)  =  S(x,  0)  =  0 
and 

P(x)=-z2l(x,0)  .  (48) 

Thus, 

ri(x,  p)  =  z21(x,  0)  -  z21(x,p)  (49) 

and  the  taper  profile  is  completely  specified  by  z2^(x, p). 

^Let  §(x,  jiu)  =  j6(x,  <u),  uu  real.  Then,  from  the  positive -slope  property 
of  reactances, 

Tg(x,w)  =  >  o  ,  (49a) 

all  uu. 
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To  exploit  the  basic  inversion  formula  (47),  it  is  at  least  necessary 
to  exhibit  a  first-order  relationship  between  the  general  coefficient  cm(p) 
and  the  segment  of  the  impulse  response  h(t)  of  the  continuous  taper  over 
the  time  interval  0<  t<  2T.  From  the  affirmative  answer  to  the  question 
posed  in  point  1  it  is  reasonable  to  infer  that  such  an  identification  exists 
and  it  is  expected  that  the  principle  developed  in  the  next  section  will  play 
an  important  role. 


THE  WEAK  CONVERGENCE  PRINCIPLE 


The  structure  of  h^(t)/R  ;  the  normalized  impulse  response  of  the 
th  ® 

n — order  discrete-cascade  approximation  to  the  continuous  taper,  is  shown 
in  Eq.  (7).  From  this  expression  it  appears  intuitively  evident  that  the 
normalized  impulse  response  h(t)/Rg  of  the  continuous  taper  must  have  the 
form 


=  6(t)  +  2u(t)  ,  (50) 

R0 

where  u(t)  is  continuous  for  t>  0  and  assumes  the  true  initial  value  u(0)  =  0.*** 

Let  i^(t)  and  i2(t)  be  any  two  smooth  input  current  drivers**  and  let 
Vj(t),  v^(t)  and  v2(t),  v^(t)  denote  the  respective  zero-state  input  voltage 
responses  at  the  terminals  of  the  continuous  taper  and  its  discrete  n—  order 
approximation.  We  make  the  consistent  physical  assumption  that  the  energy 
supplied  to  the  approximation  up  to  any  time  t  converges  as  n-»»  to  that 
supplied  to  the  continuous  taper  up  to  time  t.  Symbolically, 

t  .  ,  t 

limit  f  vj.n,(y)ik(y)dy  =  f  vk(y)ik(y)dy  ,  (51) 

n  ->  oo  o  o 

0<  t<  2t  ;  k  =  1,2 

Since  the  taper  and  the  discrete  cascade  are  linear  networks,  the 
current  i.(t)  +  i2(t)  produces  the  respective  input  voltages  v^(t)  +  v2(t)  and 
Vjn*(t)  +  vy)(t).  Thus,  (51)  also  implies  that 

limit  J  (v^+v^Hij  +  i^dy*  f  (Vj  +  v2)(ij  +  i2)dy  (52) 

n  -»oo  o  o 

and  by  subtraction  we  obtain 

limit  J  (v(1n)i2  +  v^n)i1)dy=  J  (v.i2  +  v2i  )dy  ,  (53) 

n  -»oo  o  o 

0<t<  2T  . 

Choose  i^(t)  =  6(t).  Then,  by  definition,  v^(t)  =  h^(t)  and  v^(t)  *h(t) 
so  that  (53)  yields 

In  general,  u(0+)  ^0. 

**i|(t)  and  i^(t)  are  infinitely  differentiable  for  t>0  and  all  their 
derivatives  vanish  for  t=  0. 


v^(0)  +  limit  Jh^(y)i2(y)dy  =  v_(0)  +  j  h(y)i-(y)dy  ,  (54) 

c  n-»«o  o  c  o 

0<  t<  2T  . 

Hence,  for  every  choice  of  smooth  driver  i(t), 
n  t  t 

limit  /  h  (y-2mT)i(y)dy=  J  u(y)i(y)dy,  0<t<  2T  ,  (55) 

n^»m=i  o  o 


and 

limit  v^n^(0)  =  v(0)  .  (56) 

n->oo 

In  words,  the  impulse  response  of  the  discrete  cascade  converges  weakly 
(in  the  sum  of  functional  analysis  in  Hilbert  space)  to  the  impulse  response 
of  the  continuous  taper: 

h(n)(t)-*  h(t)  ,  0<t<2T  .  (57) 

We  are  now  in  a  position  to  derive  a  time -domain  integral  equation 
which  solves  the  inverse  problem  for  the  dispersionless  taper  (§  =  0).  This 
taper  is  usually  described  by  the  telegrapher's  equation. 
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VI.  INVERSION  OF  THE  CLASSICAL  TELEGRAPH  EQUATION 
By  definition,  the  taper  is  without  dispersion  if 

*  F  *p  =  j(D 

in  uu  for  every  fixed  x;  i.e. ,  if  ?(x,  p)  is  a  constant  in  p.  However,  since 
§(x,p)  is  a  reactance  this  constant  must  equal  zero  whence, 

zn(x,p)  +  z22(x,p)  =  0  .  (59) 

It  now  follows  easily  from  the  Foster  character  of  Z(x,  p)  that 

zu(x,p)  =  z22(x,p)  ■  0  (60) 

and  that  z21(x,  p)  =  a  frequency -independent  constant.  But  then,  if  the  taper 
is  also  O.C.  transparent,  T](x,  0)  =  0  and  we  conclude  from  (48)  that  this 
constant  equals  -P(x).  To  summarize,  for  a  D.  C.  transparent  dispersionless 
taper, 

_  _ ,  A  /  Ip) 

P(x)  =  limit  (-1)  "  .  ■— r - rrr  ,  x  =  mT .  (6l) 

T*0  TAm-l(p) 

Thus,  in  the  absence  of  dispersion  the  taper  profile  is  delineated  solely 
by  its  logarithmic  taper 


and  the  discrete  approximation  is  nothing  more  than  a  cascade  of  n+1  ideal 
TEM  lines  of  the  same  1-way  delay  T, 

Consequently,  invoking  the  results  in  Problem  1,  we  find  that 

h  (t)«c  6  (t)  (63) 

where  all  the  coefficients  c^  are  real  constants,  m=  1  -*n.  Concomitantly, 
A and  Am  ^  are  independent  of  p  and  for  every  n  >  0, 

n 

=  6(t)  +  £  2cm6(t-2mT)  ,  (64) 

m=  1 


0  <t<  2(n+l) t 


Select  any  fixed  t  in  0<t<  2T  and  choose  the  integer  k  so  that 

2kt  <  t<  2(k+l)T  {65) 


where 


__  T 

t'H+T 


According  to  (55),  for  any  smooth  i(t). 


(66) 


limit  Y,  cmi(2mT)  =  J  u(y)i(y)dy  . 


(67) 


n-wo 


m=l 


But  from  the  definition  of  a  Riemann  integral. 


o 


k 

u(y)i(y)dy  =  limit  Y  2Tu(2mr)i{2mr) 
n'+0°  m=  1 


(68) 


and  a  comparison  with  (67)  immediately  leads  to  the  first-order  connection 
formula, 

cin  =  2Tu(2mT)  ,  ro  =  l-»n  .  (69) 

As  expected,  as  t  -» 0,  h^(t)  provides  an  increasingly  accurate  sampled- 
data  vefsion  of  h(t). 

Theorem  2.  Let 


6(t)  +  2u(t) 


(70) 


constitute  a  (normalized)  record  of  the  impulse  response  of  a  dispersionless 
D. C,  transparent  continuous  taper  over  the  time  interval  0<  t<  2T,  For 
fixed  x,  0<  x<  T,  let  4>(t)  denote  the  solution  of  the  integral  equation 


*M*)  +  J  u(  |t-y|)i|/(y)dy  =  6(t)  ,  (71) 

o 

0  <  t  <  2x  . 

Then, 

P(x)=  -2t(/(2x),  0  <  x  <  T  .  (72) 
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Proof.  (Largely  heuristic).  Let 


v|i(t)  =  6(t)+cp(t)  . 

Clearly,  44*)  =  t>  0,  and 

.  2x 


tp(t)+/  u(jt-yj)  cp(y)dy= -u(t)  , 


(73) 


(74) 


0  <  t  <  2x  . 

This  conventional  distribution-free  Fredholm  integral  equation  of  the  second 
kind  is  the  precise  equivalent  of  (71). 

With  T  given  by  (66)  and  the  integer  m  chosen  so  that  mT<  x<  (m+l)T, 
we  set  cp(2kT)  =  x^,  k=0-»m,  and  evaluate  the  integral  as  a  Riemann  sum  at 
the  points  t=  0,  2t,  . . ,  2m T  with  the  aid  of  (69).  The  result  is  the  linear 
system, 
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1 

C1 

C2 

• 

• 

• 

cm 

xo  +  I? 

-IS 

1 _ 

C1 

1 

C1 

• 

• 

• 

Cm-1 

X1 

0 

C2 

C1 

1 

• 

• 

• 

• 

x2 

= 

0 

• 

• 

• 

• 

• 

• 

% 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

c 

m 

Cm-1 

c  -> 
m-2 

• 

• 

• 

1 

X 

m 

—  — 

0 

(75) 


(Observe  that  the  derivation  of  (75)  has  made  use  of  the  initial  value  constraint, 
u(0) =  0. ) 


It  now  follows  easily  from  Cramer's  rule  that 

A 


-~i  =  2Tx  +1 
o 


(76) 


and 


m 


A  / 

(• « m  m 
-1)  •  ■  A  —  =  2x 

'  TA__  m 
m 


(77) 


Thus, 


■HIT 


This  algebraic  procedure  works  in  a  general  L2*setting  [4], 


m  ,  ,,m-l 
Ztx  +1  " l"1' 


A  , 
m 

TA  , 
m-1 


-2  limit  x  = -2<p(2x)  = -24j(2x)  =  limit  (-1)1 
TjO  m  t}0 


=  P(x)  ,  (79) 


Q.E.  D. 


Comment  1.  The  determination  of  P(x)  for  0<  x<  Y  has  been  reduced 
to  the  solution  of  a  1-parameter  family  of  Fredholm  integral  equations  of  the 
second  kind  (with  parameter  x).  Moreover,  because  of  the  symmetric  positive, 
definite  character  of  the  coefficient  matrix  in  (75)  and  the  sampled-data  charac¬ 
ter  of  the  cm's,  equation  (79)  can  be  made  the  basis  of  a  stable  digital  algo¬ 
rithm  of  the  Levinson  type  which  operates  in  real  time.  In  this  setting  it  is 
of  course  natural  and  (advisable)  to  exploit  the  theory  of  orthogonal  poly¬ 
nomials  on  the  unit  circle. 

As  an  illustrative  example,  consider  the  exponential  line  whose  Laplace 
transform  description  is  given  by 


=  -pR(x)I(x,  p)  , 
=-pR"1(x)V(x,p) 


where 


R(x)  =  RQe  ,  0<x<T  .  (! 

As  usual,  Rq  and  a  are  real  constants,  the  former  being  positive,  and 

v(x,t)  »V(x,p)  ,  (I 


Thus, 


i(x,t)  «I(x,p)  . 


p(x)  =  £  .  =  a  , 


Note  that  x  is  in  secs,  since  the  local  velocity  of  propgation  has  al¬ 
ready  been  normalized  to  unity. 
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a  constant 


To  calculate  the  corresponding  input  voltage  impulse  response 
v(0,  t)  =  h(t),  we  first  eliminate  I  to  obtain  the  second-order  differential 
equation  for  V, 


.2adV.p2v*o  . 

(86) 

Hence, 

V(x,  p)  =  Aeta-6)*+Beta+6):t 

(87) 

where  A(p)  and  B(p)  are  independent  of  y  and 

5(P)  =^P^+aZ  . 

(88) 

(We  choose  that  branch  of  the  square -root  for  which  6(p)>  0  for 
guarantees  that  6(p)  is  positive-real  in  Rep>  0.) 

p>  0.  This 

14 

If  the  line  is  semi -infinite,  B(p)  =  0  and  we  find  that  for  Rep>  0 

V(x,p)  =  Ae^a"6^*  , 

(89) 

(90) 

and 

z.  ^0 

in  1(0,  p)  j~2 — j 

Vp  +a  'a 

Thus,  since  h(t)  «►  Zin(p),  (91)  yields 

/  2  2 

h(t)  J P  +a  -a 

Tp-  «*  - - - -  .  Rep>0  , 

Ko  P 

(92) 

and  carrying  out  the  inversion  we  obtain  the  explicit  result 

^  =  6(t)  +  2u(t)  , 

R0 

(93) 

^The  exponential  e^a+^^x-»«  as  x  -»»  since  a+6(p)  has  positive  real- 
part  in  the  right -half  p-plane, 
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where 1 


u(t)  =  |-(l  +  /  JQ(y)dy  -  Jx(at)) 


/I  +  J.  (at)  *  \ 

=  »  (~  -f—  +  £  J2r+l<at0’ 


r=l 


(94) 

(95) 


0<  t<  2T  . 


Note  that  the  analytic  form  of  the  impulse  response  is  quite  complicated, 
even  for  this  very  simple  nonuniform  line. 

Another  excellent  and  challenging  example  is  the  trigonometric  line 
for  which 

R(x)  =  Rq/cos^x  ,  0<x<  T  .  (96) 

Clearly, 

P(x)  =  tanx  (97) 

and  the  behavior  close  to  x=  rr/2  is  most  severe. 

The  relevant  Diff.  Eq.  for  V(x,  p)  is  easily  found  to  be 


,  t  dV  2.,  n 
-  2  tan  x  •  "  P  V  =  0 


(98) 


and  its  solution  is  given  by 

—  n — 

V(x,  p)  =  (A(p)  cosh.y  p  -1  y +  B(p)sinh,yp  -ly)/cos  x,  Rep>0.  (99) 

It  is  now  a  straightforward  matter  to  derive  the  relatively  simple  expression 


^  =  6(t)  +  Ij(t),  t>  0  , 


(100) 


where  I^(t)  is  the  modified  Bessel  function  of  order  one.  Clearly, 

u(t)  =  1^/2  .  (101) 


(y)  is  the  r^  order  Bessel  function  of  the  first  kind. 
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(The  agreement  between  the  exact  solutions  (85)  and  (97)  and  those  obtained 
by  the  use  of  (72)  has  been  encouraging  but  a  further  refinement  of  the  soft* 
ware  is  still  needed. ) 

In  Part  2  of  this  report  we  hope  to  be  able  to  describe  the  complete 
solution  for  the  dispersive  case.  Our  research  up  to  this  point  appears  to 
indicate  that  the  unknown  per-unit  length  group-delay  function  ?(x, p)  intro¬ 
duces  an  additional  integral  equation.  The  final  formulation  will  probably 
involve  a  1  -parameter  family  of  two  coupled  Fredholm  integral  equations  of 
the  second  kind. 
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APPENDIX  1 


Let  the  reactance  2-port  N/  delimited  between  the  dotted  lines  in  Fig. 
be  described  by  its  2x2  scattering  matrix  S(p)  normalized  to  Rq  on  the  left 
and  on  the  right.  It  can  be  shown  by  straightforward  analysis  [3]  that 
there  exist  two  real  2 -variable  polynomials 

n 

MP.  *)  =  J  zrhr(p)  (1.1) 

r=0 


and 


g(P»  ®)  =  Yj 
r=0 


(1.2) 
16 


in  p  and  z  and  a  real  polynomial  f(p)  ^  0  in  the  single  variable  p  such  that 


S(p)  =  W(p,e"2pT)  =  W(p,  z) 


z=e‘2pT 


(1.3) 


where 


W(p,  z)  = 


1 

» 

h(p,  z) 

e"npTf{p) 

gtp.  z) 

_=‘np\(P) 

-z^tp,  z) 

(1.4) 


The  matrix  W(p,  7.)  is  2-variable  bounded-real  (b.r.).  gip,  z)  is  2- 
variable  Hurwitz  and  gQ(p)  is  Hurwitz;  i.e. ,  if  we  let 

D+(l)  =  {(p,  z):Rep>  0,  |z|<l],  (1.5) 


it  is  true  that 


and 


12-W  (p,  z)  W(p,  z)  >02 


g(p.  z)  i  o 


(1.6) 


(1.7) 


for  all  pairs  (p,  z)eD  (1)  and 


gn(p)  t  0  »  Rep>  0  . 


(1.8) 


*^The  use  of  z  as  a  symbol  for  impedance  and  as  a  complex  variable 
should  cause  no  confusion. 
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In  addition,  h(p,z),  g(p,  z)  and  £{p),  viewed  as  polynomials  in  two  abstract 
variable  p  and  z,  are  tied  together  by  the  Feldtkeller  relation 


g(P»  *)g*(p»  *>  -  h(P.  a)^(p,  z)  =  f(p)^(p)  =  v(p)^  0  . 


(1.9) 


Our  immediate  aim  is  to  exploit  the  evident  intimate  connection  between 
the  2-variable  input  reflection  coefficient 


sn(p,z>  '  gfetf 


(1.10) 


and  the  backend  input  reflection  coefficient 


s22(p,z)  =- 


a^tP.a) 

g(p.a) 


(1.11) 


This  we  accomplish  by  means  of  the  following  piece  of  physical  insight. 

Let  Zj(p,z)  and  Z^(p,  z)  denote  the  2-variable  driving -point  impedances 
seen  looking  into  the  front  and  back  ends  of  N7  with  the  opposite  ports  closed 
on  Rn  and  R^,  respectively.  Since  z  is  ultimately  to  be  identified  with  e~^PT 
which  approaches  zero  as  T  -» oo  for  all  fixed  p  in  Re  p>  0,  it  is  obvious  that 
Z2(P»  0)  is  nothing  more  than  the  backend  input  impedance  of  N7  with  all  lines 
infinitely  long.  However,  since  the  input  impedance  of  an  infinitely  long  line 
equals  its  characteristic  impedance,  this  means  th^t  Z^(p,0)  is  the  impedance 
seen  looking  into  the  backend  of  LLN  n  with  its  input,  port  closed  on 
Thus,  in  the  notation  of  (11), 


Z  (p,  0)  -  R  .  . 

s22(p,°)=  z  ip  0)  +  R  =  s22(p)  . 

c  n 

Consequently,  using  (1.1),  (1.2),  and  (1.11)  we  obtain, 


s7?(p,z)  = 


hn*(P)  +  zhn_  !»(?)+♦•«  +  »  hp»(p) 
g0(P)  +  *g  j(P)  +  •  •  •  +  *ngn\P) 


(1.12) 


(1.13) 


(nl  h„#(P> 

s22(p,  0)  =  s22  (p)  =  -  -  7pj- 


(1.14) 


a  key  result. 

By  definition, 


-24- 


(1.15) 


1+s 


-  r  1 1  _  R  g+h 

'l'R0‘  F-'s '1 1  ~  R0  *  g^TT 


and  a  direct  calculation  yields 


v 


(g-h)(g-h)*  "  (g-K)(g-h)#  "  ee# 


•where 


n  n 

e(p,z)  sg(p,z) -h(p,  z)=  Yj  zr(g r“hr)  =  Z  *r®r(P) 

r=0  r=0 

Clearly,  Z  j  (p,  0)  =  RQ  so  that 

1>0<P) 

0*an(p>0)  “i^TpI  • 


Hence, 


From  (1.9), 


and  (1.18)  forces 


since  gn(p)  ^  0.  Thus, 


hQ(p)  s  0 


8n(p)g0*(p)  "  hn(p)V(p)  £  0 


gn(p)  -  0 


eQ(p)  =  gQ(p) -hQ(p)  =  g0(p)  , 


en(p)  *  8ntP)  '  hn(P>  =  -hn(P> 


and  we  conclude  from  (1.14)  that 

(n).  .  _  en*^ 
s22(p)  -  eQ(p)  * 

Since  Z^(p,  z)/Rq  is  2-variable  p.  r.  it  admits  an  expansion 

Zl(p,z)  “  r 

- TT - >  +  2  Z  .  cp(P) 

0  r=  1 

-25- 


(1.16) 

(1.  16a) 

(1.17) 

(1.18) 

(1.19) 

(1.20) 

(1.21) 

(1.22) 

(1.23) 

(1.24) 


which  converges  for  all  z  in  |z  |<  1  for  every  fixed  p  in  Re  p>  0  and  for  all 
but  a  finite  number  of  p's  on  the  juu-axis.  These  assertions  are  easily 
validated  by  making  use  of  the  Schiir  criterion  for  functions  with  positive 
real-part  in  |z|<  1,  Namely,  for  every  r>l, 


Z(r)(p)  = 


Zc^p) 

2c2(p) 


^(p) 


2cr(p)  2cr  l(p)  .  .  .  2cl(p)  1 


is  a  rational  p.  r.  matrix  [2], 
According  to  (1, 16), 


00  00 

(1  +  Z  zrcr(P)  +  Z  z‘rcr#(p))(e0(p>+Zel(p)+**'+znen(p))= 


and  by  equating  coefficients  of  z  on  both  sides,  r  =  0-*n,  we  derive  the 
linear  matrix  equation. 


where 


(1.29) 


From  Cramer's  rule. 


and 


(n\  -  v(p)  ^n-1^ 

0(p) '  ■  “S^pT 

,  v  ,  ,4n  v(p)  ^ n  ^ 


(1.30) 


Thus,  invoking  (1.23), 


In),  ,  %»<!’> 

322<P»  =  "05T 


,.»  e0»(p>  V*(p) 


(1.31) 


Now  if  we  go  back  to  Eq.  (12)  and  substitute  (1.31)  we  find  the  expression 

n  e0(P>  A„'(P> 

.n(p)=(-D  •  d„(p)  •  ^51  •  U.32) 

for  the  junction  reflection  coefficient  of  LLN  n. 

In  view  of  (1.29),  e^( p)  must  be  determined  from  the  equation 

An-l(p) 

eo(p)eQi!t(p) =  v(p)  •  ^  •  (1.33) 


Since  e^(p)  =  g^p)  is  Hurwitz,  it  coincides  with  the  Wiener  factor  of  the 
right-hand  side  of  (1.33).  In  general,  therefore,  Cq(p)  depends  not  only  on 
the  coefficients  Cj(p),  c2(p),  .. .,  cR(p)  but  also  on  that  part  of  v(p)  not 
cancelled  by  the  even  (nonnegative)  rational  function  Ar  ^(p)/An(p). 

To  examine  the  matter  more  closely,  let 


Vl<p>  an(p)an*lp) 

"a^pT  =  Sn(p)Bn.(p) 


(1.34) 


where  <*n(p)  and  ?n(p)  are  two  relatively  prime  Hurwitz  polynomials.  Evidently, 
^n^n#  c^vl<les  v  because  the  product  ege^  is  polynomial.  More  strongly  [3], 
if  all  LLN's  are  all-pass  free  on  their  output  sides,  v(p)  =  ©n<p) Pn*(p)  and 


en(p)  =  a„(P>  ; 


(1.35) 


i.e.,  the  Hurwitz  polynomial  egtp)  is  formed  with  all  the  zeros  of  the  quotient 
An_i(p)/An(p)  in  Re  p<  0. 

Although  we  omit  the  proof,  it  is  easily  shown  by  induction  [3]  that 
without  restriction. 


eQ<P> 


(-!)“•  n  d'^p) 
r=  1 


This  granted,  (1.32)  then  yields 


«n(p) 


An'*p> 

An-l(pl 


n- 1 


n  d'^p) 

r=l 


(1.36) 


(1.37) 


Of  course,  the  above  derivation  applies  for  all  m=  l-*n  and  we  recover  the 
diagnostic  formula  (13)  of  theorem  1. 
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